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7.1 Dark matter particles: relics from the pre-BBN era 

A general class of candidates for non-baryonic cold dark matter are weakly 
interacting massive particles (WIMPs). The interest in WIMPs as dark 
matter candidates stems from the fact that WIMPs in chemical equilibrium 
in the early universe naturally have the right abundance to be cold dark 
matter. Moreover, the same interactions that give the right WIMP density 
make the detection of WIMPs possible. The latter aspect is important as it 
provides a means to test the WIMP hypothesis. 

The argument showing that WIMPs are good dark matter candidates is 
old @, I, I, I, Ej. The density per comoving volume of non relativistic 
particles in equilibrium in the early Universe decreases exponentially with 
decreasing temperature, due to the Boltzmann factor, until the reactions 
which change the particle number become ineffective. At this point, when 
the annihilation rate becomes smaller than the Hubble expansion rate, the 
WIMP number per comoving volume becomes constant. This moment of 
chemical decoupling or freeze-out happens later, i.e. at smaller WIMP den- 
sities, for larger WIMP annihilation cross section cr a nn- If there is no subse- 
quent change of entropy in matter plus radiation, the present relic density 
of WIMPs is approximately 

ntf * 3xl0 ' 27cm3/8 . (7.1) 

\0"annf) 

For weak cross sections this gives the right order of magnitude of the DM 
density (and a temperature Tj . ~ to/20 at freeze-out for a WIMP of mass 
to). This is a ballpark argument. A more precise derivation will be presented 
in Section [721 

It is important to realize that the determination of the WIMP relic density 
depends on the history of the Universe before Big Bang Nucleosynthesis 
(BBN), an epoch from which we have no data. BBN (200 s after the Big 
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Bang, T ~ 0.8 MeV) is the earliest episode from which we have a trace, 
namely the abundance of light elements D, 4 He and 7 Li. The next observable 
in time is the Cosmic Microwave Background radiation (produced 3.8 x 10 4 
yr after the Big Bang, at T ~ eV) and the next one is the Large Scale 
Structure of the Universe. WIMPs have their number fixed at 7y o . ~ m/20, 
thus WIMPs with m Z 100 MeV would freeze out at T Z 4 MeV and would 
thus be the earliest remnants. If discovered, they would for the first time 
give information on the pre-BBN phase of the Universe. 

As things stand now, to compute the WIMP relic density we must make 
assumptions about the pre-BBN epoch. The standard computation of the 
relic density relies on the assumptions that the entropy of matter and radi- 
ation was conserved, that WIMPs were produced thermally, i.e. via interac- 
tions with the particles in the plasma, that they decoupled while the Uni- 
verse expansion was dominated by radiation, and that they were in kinetic 
and chemical equilibrium before they decoupled. These are just assump- 
tions, which do not hold in all cosmological models. In particular, in order 
for BBN and all the subsequent history of the Universe to proceed as usual, 
it is enough that the earliest and highest temperature during the radiation 
dominated period, the so called reheating temperature Trh, is larger than 



4 MeV [55(. At temperatures higher than 4 MeV, when the WIMP freeze 
out is expected to occur, the content and expansion history of the Universe 
may differ from the standard assumptions. In non-standard cosmological 
models, the WIMP relic abundance may be higher or lower than the stan- 
dard abundance. The density may be decreased by reducing the rate of 
thermal production (through a low Trjj < Tf ) or by producing radiation 
after freeze-out (entropy dilution). The density may also be increased by 
creating WIMPs from decays of particles or extended objects (non-thermal 
production) or by increasing the expansion rate of the Universe at the time 
of freeze-out. 

Non-thermal production mechanisms may also be at work within standard 
cosmological scenarios. For example, WIMPs may be produced in the out 
of equilibrium decay of other particles whose density may be fixed by ther- 
mal processes. A particular type of heavy WIMPs, WIMPZILLAS, could 
be formed during the reheating phase at the end of an inflationary period 
through gravitational interactions. Another production mechanism involves 
quantum- mechanical oscillations, for example a sterile neutrino may be pro- 
duced in the early universe by the oscillation of active (interacting) neutrinos 
into sterile neutrinos (for the latter, see the chapter of M. Shaposhnikov in 
this book). 
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In the rest of the Chapter we review the standard production mechanism 
and some of the non-standard scenarios. 

7.2 Thermal production in the standard cosmology 

7.2.1 The standard production mechanism 

In the standard scenario, it is assumed that in the early universe WIMPs 
were produced in collisions between particles of the thermal plasma during 
the radiation dominated era. Important reactions were the production and 
annihilation of WIMP pairs in particle-antiparticle collisions, such as 

XX <* e+e" , n + n~,qq, W + W~, ZZ, HH,.... (7.2) 

At temperatures much higher than the WIMP mass, T 3> m x , the col- 
liding particle-antiparticle pairs in the plasma had enough energy to create 
WIMP pairs efficiently Also, the inverse reactions that convert WIMPs into 
standard model particles (annihilation) were initially in equilibrium with the 
WIMP-producing processes. Their common rate was given by 

Tann = (o"annV)ra e q, (7.3) 

where a ann is the WIMP annihilation cross section, v is the relative velocity 
of the annihilating WIMPs, n eq is the WIMP number density in chemi- 
cal equilibrium, and the angle brackets denote an average over the WIMP 
thermal distribution. 

As the universe expanded, the temperature of the plasma became smaller 
than the WIMP mass. While annihilation and production reactions re- 
mained in equilibrium, the number of WIMPs produced decreased expo- 
nentially as e~ mx l T (the Boltzmann factor), since only particle-antiparticle 
collisions with kinetic energy in the tail of the Boltzmann distribution had 
enough energy to produce WIMP pairs. At the same time, the expansion 
of the universe decreased the number density of particles n, and with it the 
production and annihilation rates, which are proportional to n. When the 
WIMP annihilation rate r ann became smaller than the expansion rate of the 
universe H , or equivalently the mean free path for WIMP-producing colli- 
sions became longer than the Hubble radius, production of WIMPs ceased 
(chemical decoupling). After this, the number of WIMPs in a comoving 
volume remained approximately constant (or in other words, their number 
density decreased inversely with volume). 

In many of the current theories, WIMPs are their own antiparticles. For 
this kind of WIMPs (e.g. neutralinos and Majorana neutrinos), the WIMP 
density is necessarily equal to the antiWIMP density. In the following we 
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restrict our discussion to this case. We refer the reader interested in cos- 
mological WIMP-antiWIMP asymmetries, as might apply for example to a 



Dirac neutrino, to 10| 



The current density of WIMPs can be computed by means of the rate 
equation for the WIMP number density n and the law of entropy conserva- 
tion: 

^ = -3Hn - (cr ann f ) (n 2 - n 2 q ) , (7.4) 

Here t is time, s is the entropy density, H is the Hubble parameter, and as 
before n eq and {cr ann v} are the WIMP equilibrium number density and the 
thermally averaged total annihilation cross section. The first and the second 
term on the right hand side of Eq. 17.41 take into account the expansion of the 
Universe and the change in number density due to annihilations and inverse 
annihilations, respectively. 

It is customary (see e.g. @, E3, S E3]) to combine Eqs. Q7ID and ([73]) 
into a single one for Y = n/s, and to use x = m/T, with T the photon 
temperature, as the independent variable instead of time. This gives: 

f = 55|W <"■-*>■ <™> 

Here and in the rest of the Chapter we will simply write (crv) for {a ann v) 
when no ambiguity can arise. 

According to the Friedman equation, the Hubble parameter is determined 
by the mass-energy density p as 



where M P = 1.22 x 10 19 GeV is the Planck mass. The energy and entropy 
densities are related to the photon temperature by the equations 

7T 2 9-7T 2 

P=- 9cS (T)T\ s = —h cS (T)T*, (7.8) 

where g e s(T) and h e g(T) are effective degrees of freedom for the energy 
density and entropy density respectively. Recent computations of g e s(T) 
and h c g(T) that include QCD effects can be found in Ref. (5?]]. If the 
degrees of freedom parameter gl^ 2 is defined as 

1/2 _ hcs_ ( ^ _j_ 1 T dh eS \ 

1/2 

9 c s 
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then Eq. (j7.6[) can be written in the following way, 

This single equation is then numerically solved with the initial condition 
Y = Y eq at x ~ 1 to obtain the present WIMP abundance Yq. From it, the 
WIMP relic density can be computed as 

n x h 2 = ^ = m ^° h2 = 2.755 x 10 8 Y m x /GeV , (7.11) 

Pc Pc 

where and sq are the present critical density and entropy density respec- 
tively. In obtaining the numerical value in Eq. (|7.1ip we used To = 2.726 K 
for the present background radiation temperature and /i e ff(^o) = 3.91 cor- 
responding to photons and three species of neutrinos. 

The numerical solution of Eq. (|7.10p . see Fig. 17. II for an illustration, shows 
that at high temperatures Y closely tracks its equilibrium value Y eq . In fact, 
the interaction rate of WIMPs is strong enough to keep them in thermal and 
chemical equilibrium with the plasma. But as the temperature decreases, 




Fig. 7.1. Typical evolution of the WIMP number density in the early universe 
during the epoch of WIMP chemical decoupling (freeze-out). 
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Y eq becomes exponentially suppressed and Y is no longer able to track its 
equilibrium value. At the freeze out temperature (Tf. Q .), when the WIMP 
annihilation rate becomes of the order of the Hubble expansion rate, WIMP 
production becomes negligible and the WIMP abundance per comoving vol- 
ume reaches its final value. In the standard cosmological scenario, the WIMP 
freeze out temperature is about Tf. , ~ m x /20, which corresponds to a typ- 
ical WIMP speed at freeze-out of v to . = (3T f . ./2m x ) 1 / 2 ~ 0.27c. 

An important property that Fig. [7j] illustrates is that smaller annihilation 
cross sections lead to larger relic densities ("The weakest wins.") This can 
be understood from the fact that WIMPs with stronger interactions remain 
in chemical equilibrium for a longer time, and hence decouple when the 
universe is colder, wherefore their density is further suppressed by a smaller 
Boltzmann factor. This leads to the inverse relation between Q x h 2 and 
(a ann v) in Eq. 

From this discussion follows that the freeze out temperature plays a promi- 
nent role in determining the WIMP relic density. In general, however, the 
freeze out temperature depends not only on the mass and interactions of 
the WIMP but also, through the Hubble parameter, on the content of the 
Universe. Some examples of how modifications of the Hubble parameter 
affect the WIMP density are discussed in Section 17.41 below. 



7.2.2 Annihilations and coannihilations 

A pedagogical example of the dependence of the relic density on the WIMP 
mass is provided by a thermally-produced fourth-generation Dirac neutrino 
v with Standard Model interactions and no lepton asymmetry, although it is 
excluded as a cold dark matter candidate by a combination of LEP and direct 
detection limits 17, 6sJ 68j |. Figure [?T2l summarizes its relic density Q, v h? as a 



function of its mass m u . The narrow band between the horizontal lines is the 
current cosmological measurement of the cold dark matter density f^cdm^ 2 = 
0.1131±0.0034 [75]. Neutrinos with $7 X > Sl c dm are said to be overadundant, 
those with $7 X < $7 c dm are called under abundant. A neutrino lighter than 
~ 1 MeV freezes out while relativistic. If it is so light to be still relativistic 
today (m„ < 0.1 meV), its relic density is p v = 7ir 2 T^ /120. If it was massive 
enough to have become non-relativistic after freeze out, its relic density 
is determined by its equilibrium number density as p u = m u 3C(3)T^ /2ir 2 . 
Here T v = (S/ll) 1 / 3 ^, where T 7 = 2.725 ± 0.002K is the cosmic microwave 
background temperature. A neutrino heavier than ~ 1 MeV freezes out 
while non-relativistic. Its relic density is determined by its annihilation cross 
section, as in Eq. (|7.ip . The shape of the relic density curve above ~ 1 MeV 
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Fig. 7.2. Relic density Sl v h 2 of a thermal Dirac neutrino with standard-model inter- 
actions as a function of the neutrino mass m v (solid line) . The very close horizontal 
dashed lines enclose the current la band for the cold dark matter density (75j. 



in Figure 17.21 is a reflection of the behavior of the annihilation cross section 
into lepton-antilepton and quark-antiquark pairs //: the Z-boson resonance 
at m„ ~ mz/2 gives rise to the characteristic V shape in the Qh 2 curve. 
Above m u ~ 100 GeV, new annihilation channels into Z- or W-boson pairs 
open up (thresholds at m n u = mw and m u = mz, respectively). When the 
new channel annihilation cross sections dominate the relic density decreases. 
Soon, however, the perturbative expansion of the cross section in powers of 
the (Yukawa) coupling constant becomes untrustworthy (the question mark 



in Figure 173]) . A general unitarity argument 16] limits the relic density to 
the dashed curve on the right, 



The relic density of other WIMP candidates exhibits features similar to 
that of the Dirac neutrino just discussed. In general, because of the presence 
of resonances and thresholds in the annihilation cross section, one should 
not rely on a Taylor expansion of cr ann ?; in powers of v, because it would 
lead to unphysical negative cross sections. Let us remark that resonant 
and threshold annihilation, including the coannihilation thresholds discussed 
below, are ubiquitous for neutralino dark matter (see the of chapter of J. Ellis 
in this book). In the non-rGl&tivistic limit v — )■ 0, th.6 product o" a nn^ 

tends to 

a constant, because of the exoenergetic character of the annihilation process 
that makes the annihilation cross section <r ann diverge as 1/v as v — >■ 0. One 
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can safely Taylor expand a ann v in powers of v 2 if a ann v varies slowly with v, 

CannV = a + bv 2 H . Then the thermal average is (o" ann u) = a + bM- + • • • . 

Close to resonances and thresholds, however, <7 ann varies rapidly with v 
and more sophisticated procedures (described in the following) should be 
used 



2l|,|20(. 



State-of-the-art calculations of WIMP relic densities strive to achieve a 
precision comparable to the observational one, which is currently around 
a few percent. Since the speed of the WIMPs at freeze-out is about c/3, 
relativistic corrections must be included. Fully relativistic formulas for any 
cross section, with or without resonances, thresholds, and coannihilations, 



were obtained in |20l. 1271]. In the simplest case without coannihilations, one 
has 

_ I^lw^lM^IIl, (7 . 13) 

where W xx (s) is the XX annihilation rate per unit volume and unit time (a 
relativistic invariant), s = 4(m x +p 2 ) is the center-of-mass energy squared, 
and Ki(x), K2( x ) are modified Bessel functions. The Lorentz-invariant 
annihilation rate Wij(p) for the collision of two particles of 4- momenta pi 
and pj is related to the annihilation cross section a™ through 

W ij (s) = a ij F ij , (7.14) 

where 



F- ■ 



(Pi - Pj) 2 - mfm 2 (7.15) 



is the Lorentz-invariant flux factor. 

Coannihilations are an essential ingredient in the calculation of the WIMP 
relic density. They are processes that deplete the number of WIMPs through 
a chain of reactions, and occur when another particle is close in mass to the 
dark matter WIMP (mass difference Am ~ temperature T). In this case, 
scattering of the WIMP off a particle in the thermal 'soup' can convert 
the WIMP into the slightly heavier particle, since the energy barrier that 
would otherwise prevent it (i.e. the mass difference) is easily overcome. The 
particle participating in the coannihilation may then decay and/or react 
with other particles and eventually effect the disappearance of WIMPs. We 
give two examples in the context of the Minimal Supersymmetric Standard 
Model. Neutralino coannihilation with charginos x* 1 may proceed, for in- 
stance, through 

X?e-^x 2 ^ e , X 2 ^X2<lu, X2Xi^W+W~. (7.16) 
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Fig. 7.3. The effective invariant annihilation rate W c s as a function of p e g for a 
particular supersymmetric model examined in [27j . The final state threshold for 
annihilation into W + W~ and the coannihilation thresholds appearing in Eq. (|7.19l) 
are indicated. The X2X2 coannihilation threshold is too small to be seen. 



Neutralino coannihilation with tau sleptons f may instead involve the pro- 
cesses 



T7, 



rXx 



(7.17) 



Coannihilations were first included in the study of near-degenerate heavy 



neutrinos in |8J and were brought to general attention in [211 ] . The rela- 
tivistic treatment was formulated in [271 ] ■ Under the conditions described 
below, which are reasonable during WIMP freeze-out, one replaces (cr ann v} 
in Eq. ([73} with 



/ °° d PeS p 2 eS W eS (s) K x {^fs/T) 



4- 



1 2 ' 



(7.18) 



where s = 4p e g + 4m?,, Oj is the number of internal degree of freedom 
(statistical weight factor) for the i-th coannihilating particle, and 

Fij fli'lj 



W eS (s) = Y, 



F xx 9 



2 Wijis). 



(7.19) 



The sums extend over all the N coannihilating particles, including the x-> 
and mi = m x , g\ = g x . The assumptions underlying Eq. (|7.18p are: (1) 
all coannihilating particles decay into the lightest one, which is stable, and 



10 



G. Gelmini, P. Gondolo 



their decay rate is much faster than the expansion rate of the universe - so 
the final WIMP abundance is simply described by the sum of the density of 
all coannihilating particles; (2) the scattering cross sections of coannihilat- 
ing particles off the thermal background are of the same order of magnitude 
as their annihilation cross sections - since the relativistic background par- 
ticle density is much larger than each of the non-relativistic coannihilating 
particle densities, the scattering rate is much faster and the momentum 
distributions of the coannihilating particles remain in thermal equilibrium; 
(3) all coannihilating particle are semi-relativistic, so the Fermi-Dirac and 
Bose-Einstein thermal distributions can be replaced by Maxwell-Boltzmann 
distribution /j = e~ E ^ T . 

An important aspect of the effective annihilation rate in Eq. (|7.18p is that 
coannihilations appear as thresholds at a value of y/s equal to the sum of 
the masses of the coannihilating particles. As an example of this, Fig. 17.31 
taken from Ref. 27]], shows that coannihilation thresholds and regular final 
state thresholds appear on the same footing in the invariant annihilation 
rate W c g. 

Computations of WIMP relic densities can become quite involved, es- 
pecially in the presence of coannihilations. There exist publicly available 
software 5J, [65[ that can handle these calculations for generic WIMPs (see 
the chapter of F. Boudjemain in this book). 



7.3 Non-thermal production in the standard cosmology 
7.3.1 Gravitational mechanisms 



WIMPZILLAs [3J, |2& |3J, |23, |3Ji |3J illustrate a fascinating idea for gener- 
ating matter in the expanding universe: the gravitational creation of matter 
in an accelerated expansion. This mechanism is analogous to the produc- 
tion of Hawking radiation around a black hole, and of Unruh radiation in 
an accelerated reference frame. 

WIMPZILLAs are very massive relics from the Big Bang: they can be 
the dark matter in the universe if their mass is ~ 10 13 GeV. They might be 
produced at the end of inflation through a variety of possible mechanisms: 
gravitationally, during preheating, during reheating, in bubble collisions. It 
is possible that their relic abundance does not depend on their interaction 
strength but only on their mass, giving great freedom in their phenomenol- 
ogy. To be the dark matter today, they are assumed to be stable or to have 
a lifetime of the order of the age of the universe. In the latter case, their 
decay products may give rise to the highest energy cosmic rays. 

Gravitational production of particles is an important phenomenon that 
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T= H/2n 




Fig. 7.4. Relic density of gravitationally-produced WIMPZILLAs as a function of 
their mass Mx- Hi is the Hubble parameter at the end of inflation, T r h is the 
reheating temperature, and M p \ 3 x 10 19 GeV is the Planck mass. The dashed 
and solid lines correspond to inflationary models that smoothly end into a radiation 
or matter dominated epoch, respectively. The dotted line is a thermal distribution 
at the Gibbons-Hawking temperature T = Hj/2tt. Outside the 'thermalization 
region' WIMPZILLAs cannot reach thermal equilibrium. (Figure froml32h 



is worth describing here. Consider a scalar field (particle) X of mass Mx 
in the expanding universe. Let r\ be the conformal time and a{rf) the time 
dependence of the expansion scale factor. Assume for simplicity that the 
universe is flat. The scalar field X can be expanded in spatial Fourier modes 
as 

(27r)3/2Q(7?) [a k h k (v)e ik - x + aihUr,)e- ik - x ] . (7.20) 

Here a/% and at are creation and annihilation operators, and hk(rj) are mode 
functions that satisfy (a) the normalization condition h^h'^ — h' k li^ = i (a 
prime indicates a derivative with respect to conformal time), and (b) the 
mode equation 

hUv)+^k(v)hk(v) =0, (7.21) 



where 



uHn) k z + Mia 2 + (6£ - 1) — . (7.22) 

a 



The parameter £ is £ = for a minimally-coupled field and £ = i for 
a conformally-coupled field. The mode equation, Eq. (|7.2ip . is formally 
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the same as the equation of motion of a harmonic oscillator with time- 
varying frequency i^k(v)- For a given complete set of positive- frequency 
solutions hk{f]), the vacuum 10^) of the field X, i.e. the state with no X 
particles, is defined as the state that satisfies a^O^) = for all k. Since 
Eq. (|7.2ip is a second order equation and the frequency depends on time, the 
normalization condition is in general not sufficient to specify the positive- 
frequency modes uniquely, contrary to the case of constant frequency ojq 
for which h^ijf) = e~ tulori / (2ujq) 1 ^ 2 . Different boundary conditions for the 
solutions h k (77) define in general different creation and annihilation operators 
aj, and at, and thus in general different vacua!]] For example, solutions which 
satisfy the condition of having only positive- frequencies in the distant past, 

h(rj) ~ e'^k v for v (7.23) 

contain both positive and negative frequencies in the distant future, 

h(r]) ~ a k e~^kV + (3 k e +iuJ *v for v +oc . (7.24) 

Here uj^ = lim^ioo Wk(v)- As a consequence, an initial vacuum state is 
no longer a vacuum state at later times, i.e. particles are created. The 
number density of particles is given in terms of the Bogolubov coefficient f3k 
in Eq. (TMl by 



These ideas have been applied to gravitational particle creation at the end 



of inflation by 32J and 29|]. Particles with masses Mx of the order of the 
Hubble parameter at the end of inflation, Hj ~ 10~ 6 Mpi ~ 10 13 GeV, may 
have been created with a density which today may be comparable to the crit- 
ical density. Figure [731 shows the relic density Qxh 2 of these WIMPZILLAs 
as a function of their mass Mx in units of Hj. Curves are shown for infla- 
tion models that have a smooth transition to a radiation dominated epoch 
(dashed line) and a matter dominated epoch (solid line). The third curve 
(dotted line) shows the thermal particle density at temperature T = Hj/2ir. 
Also shown in the figure is the region where WIMPZILLAs are thermal 
relics. It is clear that it is possible for dark matter to be in the form of 
heavy WIMPZILLAs generated gravitationally at the end of inflation. 



f The precise definition of a vacuum in a curved space-time is still subject to some ambiguities. 
We refer the interested reader to [f| [T3 . 1771 . |78| and to the discussion in [46(1 and references 
therein. 
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7.3.2 Decays 

Dark matter may be produced in the decay of other particles. If the DM 
particles are non-interacting when the decay occurs, they inherit (except for 
some entropy dilution factor) the density of the parent particle P 

n DM h 2 ~ ^MOp h 2 . (7.26) 
mp 

This is the case of superWIMPs (see the chapter J. Feng in this book), 
extremely weakly interacting particles produced in the late decays of WIMPs 
(e.g. axinos or gravitinos from the decay of neutralinos or sleptons) which 
practically only interact gravitationally and cannot be directly detected. In 
some models the superWIMP may produce WIMPs through its decay. This 
is the case, for example, of gravitinos producing Winos (which otherwise 
would have a very low thermal relic density) with the right DM abundance 



through their decay 12], |34j . 



7.4 Thermal and non-thermal production in non-standard 

cosmologies 

The relic density (and also the velocity distribution before structure forma- 
tion) of WIMPs and other DM candidates such as heavy sterile neutrinos 
and axions, depends on the characteristics of the Universe (expansion rate, 
composition, etc.) immediately before BBN, i.e. at temperatures T > 4 



MeV [55J . The standard computation of relic densities relies on the assump- 
tion that radiation domination began before the main epoch of production 
of the relics and that the entropy of matter and radiation has been con- 
served during and after this epoch. Any modification of these assumptions 
would lead to different relic density values. Any extra contribution to the 
energy density of the Universe would increase the Hubble expansion rate 
H and lead to larger relic densities (since the decreasing interaction rate T 
becomes smaller than H earlier, when densities are larger). This can happen 
in the Brans-Dicke-Jordan 18] cosmological model, models with anisotropic 
expansion 0, Q3, H^], scalar-tensor 30, 33, 51, 71] or kination 5(], II] mod- 
els and other models [HI, 64, 70] In some scalar-tensor models H may be 



decreased, leading to smaller relic densities [71]. These models alter the 
thermal evolution of the Universe without an extra entropy production. 

Not only the value of H but the dependence of the temperature T on 
the scale factor of the Universe could be different, if entropy in matter 
and radiation is produced. This is the case if a scalar field <f) oscillating 
around its true minimum while decaying is the dominant component of 
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the Universe just before BBN. This field may be an inflaton or another 
late decaying field, such as a modulus in supersymmetric models. Models 
of this type include some with moduli fields, either the Polonyi field 



251 ] or others [36[ ) or an Affleck-Dine field and Q-ball decay [4J, [52| , and 



24, 



thermal inflation [26]. Moduli fields correspond to flat directions in the 
supersymmetric potential, which are lifted by the same mechanisms that give 
mass to the supersymmetric particles of the order of a few to 10's of TeV, 
and they usually have interactions of gravitational strength. The decays 
of the (j) field finally reheat the Universe to a low reheating temperature 
Trh, which could be not much larger than 5 MeV. In these low temperature 
reheating (LTR) models there can be direct production of DM relics in the 
decay of <fi which increase the relic density, and there is entropy generation, 
through the decay of 4> into radiation, which suppresses the relic abundance. 



Thus, in non-standard cosmological scenarios, the relic density of WIMPs 
Q x may be larger or smaller than in standard cosmologies f2 s td- The density 
may be decreased by reducing the rate of thermal production (through a low 
Trh < Tf M _ ) , by reducing the expansion rate of the Universe at freeze-out or 
by producing radiation after freeze-out (entropy dilution). The density may 
be increased by creating WIMPs from particle (or extended objects) decay 
(non-thermal production) or by increasing the expansion rate of the Universe 
at freeze-out. Usually these scenarios contain additional parameters that can 
be adjusted to modify the WIMP relic density. However these are due to 
physics that does not manifest itself in accelerator or detection experiments. 



Let us comment that not only the relic density of WIMPs but their char- 
acteristic speed before structure formation in the Universe can differ in 
standard and non-standard pre-BBN cosmological models. If kinetic de- 
coupling (the moment when the exchange of momentum between WIMPs 
and radiation ceases to be effective) happens during the reheating phase of 
LTR models, WIMPs can have much smaller characteristic speeds, i.e. be 
much "colder" (7^], with free-streaming lengths several orders of magni- 
tude smaller than in the standard scenario. Much smaller DM structures 
could thus be formed, a fraction of which may persist as minihaloes within 
our galaxy and be detected in indirect DM searches. The signature would 
be a much larger boost factor of the annihilation signal than expected in 
standard cosmologies for a particular WIMP candidate. WIMPs may in- 
stead be much "hotter" than in standard cosmologies too, they may even be 
warm DM instead of cold, which would leave an imprint on the large scale 
structure spectrum [4l|, |4£J, |62j . 
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7.4-1 Low temperature reheating (LTR) models 

Let us consider a late decaying scalar field (j> of mass and decay width 
which dominates the energy density of the Universe while oscillating 
about the minimum of its potential and decays reheating the Universe to 
a low reheating temperature Irh, with 5 MeV < Trh ^ Tf.o. for Trh, 
so BBN is not affected. The usual choice for the parameter Trh is the 
temperature the Universe would attain under the assumption that the (ft 
decay and subsequent thermalization are instantaneous, 



r„ 



Hd. 



ecay 




90 M M P 



(7.27) 



Here, is the decay width of the 4> field, ~ m^/A^ ff . If <fi has non- 
suppressed gravitational couplings, as is usually the case for moduli fields, 
the effective energy scale A g ~ Mp (but A e g could be smaller [45|]). Thus, 
with g* ~ 10, 



T RH 10 MeV 



100 TeV 



3/2 (M P 



A, 



off 



(7.28) 



Numerical calculations in which the approximation of instantaneous decay 
is not made show that the parameter Trh provides a good estimate of the 
first temperature of the radiations dominated epoch (see Fig. 17. 5|) . 

Both thermal and non thermal pro duction mechanisms in LTR modesl 

H, H 0, SB, H Si, 0, EI Si, S EL E" 



have been discussed 



fid 

Ififil], mostly in supersymmetric models where the WIMP is the neutralino. 
The decay of (j> into radiation increases the entropy, diluting the WIMP 
number density. The decay of (ft into WIMPs increases the WIMP number 
density. In supersymmetric models cp decays into supersymmetric particles, 
which eventually decay into the lightest such particles (the LSP, typically a 
neutralino). Call b the net number of WIMPs produced on average per <fi 
decay, which is a highly model dependent parameter 

A combination of Trh and the ratio b/m^ can bring the relic WIMP 
density to the desired value O ct } m [fi3]- The equations which describe the 
evolution of the Universe depend only on the combination b/m^ and not on 
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b and separately. They are 

p = -3H(p + p) + T (j> p^ (7.29) 



h x = -3Hn x - (av) (n x - n^ eq ) + — T^p^ (7.30) 
Ptj, = SHpcf, - T^p^ (7.31) 

1,1 = 3^<" + ^' < 7 - 32 ' 

In Eqs. (j7.29H7.32]) . a dot indicates a time derivative, p^ is the energy density 
in the <f> field, which is assumed to behave like non-relativistic matter; p 
and p are the total energy density and pressure of matter and radiation at 
temperature T; n x is the number density of WIMPs (which are assumed 
to be in kinetic but not necessarily chemical equilibrium) and n x>e q is its 
value in chemical equilibrium; finally, H = a/a is the Hubble parameter, 
with a the scale factor. The first principle of thermodynamics in the form 
d(pa 3 ) + d{p ( j ) a 3 ) + pda 3 = Td(sa 3 ) can be used to rewrite Eq. (|7.29p as 

s = -3Hs + ^. (7.33) 

where s = (p + p — m x n x )/T is the entropy density of the matter and 
radiation. For p^ — > these equation reduce to the standard scenario. 
During the </>-oscillation-dominated epoch, H oc T 4 22|. This can be seen 



using Eq. (|7.29|) while the matter content is negligible. In Eq. (|7.29p with 
p = p/S) substitute p ~ T 4 and p^ ~ MpH 2 . Then use H ~ t , write 
T oc t a , where a is a constant, match the powers of t in all terms, and 
determine that a = —(1/4). Hence, H oc t^ 1 oc T 4 (and p^ oc H 2 oc T 8 ). 
Since H equals 1r H /M p at T = T RH , it is H ~ T 4 /(T| H M P ). 

The initial conditions are specified through the value Hj of the Hubble 
parameter at the beginning of the (^-oscillations dominated epoch. This 
amounts to giving the initial energy density p^j in the <f> field at the be- 
ginning of the reheating phase, or equivalently the maximum temperature 
of the radiation Tmax- Indeed, one has Hi ~ p\j/Mp ~ ^max/ '(^rh-^p)- 
The latter relation can be derived from p^ ~ T 8 /T^ H and the consideration 
that the maximum energy in the radiation equals the initial (maximum) 
energy p^j. As the <fi begins to decay, the temperature of the radiation bath 
rises sharply to Tmax [3J|, decreases slowly as function of the scale factor 
a during the ^-oscillating dominated phase, as T ~ a -3 / 8 until it reaches 
7rh> when the radiation dominated phase starts and T ~ a -1 . 

Fig. 17.5b , shows how the WIMP density £l x h 2 depends on Trh for illustra- 
tive values of the parameter r\ = ^(lOOTeV/m^), both for WIMPs which are 
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Fig. 7.5. a. (top) WIMP density ft x h 2 as function of the reheating temperature Trh 
for illustrative values of the ratio r\ = ^(lOOTeV/m^) [63| . b. (middle) Evolution of 
the neutralino x abundance for different values of Trh and rj = in an mSUGRA 
model with M 1/2 = m Q = 600GeV, A = 0, tan/3 = 10, n > 0, m x = 246GeV and 
standard relic density Q s tdh 2 — 3.6 [6l|. The short vertical lines indicate T^h- [6l| . 
c. (bottom) Same as b. but for Trh = 1 GeV and several values of r/. 
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underdense and for WIMPs that are overdense in usual cosmologies. The 
behavior of the relic density as a function of Trh is easy to understand phys- 
ically. The usual thermal production scenario occurs for Trh > T^ . For 
Trh < 7f ., there are four different ways in which the density Qh 2 depends 
on Trh- There are four cases (6^]: (1) Thermal production without chem- 



ical equilibrium, for which Q x ~ T^ H 31]. (2) Thermal production with 
chemical equilibrium, in which case the WIMP freezes out while the uni- 
verse is dominated by the (p field. Its freeze-out density is larger than usual, 
but it is diluted by the entropy produced in <j) decays (Fig. 17.5b ). In this 
case f2 x oc 7r H . (3) Non-thermal production without chemical equilibrium, 
where £l x oc ??Trh (independently of any assumption on neutralino kinetic 
equilibrium) (Fig. 17.5b ) . (4) Non-thermal production with chemical equilib- 
rium, where Q oc (Fig. 17.5b ). For the validity of the annihilation term 
in Eq J7.29I one needs to assume that WIMPs enter into kinetic equilibrium 
before production ceases. In any event the solutions just presented should 
remain qualitatively valid because kinetic equilibrium affects only solutions 
which interpolate in Trh between two correct solutions, namely the solution 
of the standard cosmology at high Trh for which WIMPs are initially in 
kinetic equilibrium, and the WIMP production purely through the scalar 
field decay (case 3), for which kinetic equilibrium is irrelevant. 

For all overdense (f^td > ^cdm) WIMPs, given one value of r] < 10 
(100GeV/m x ) there is only one value of Trh for which Q x = Cl c dm- The 
exception is a severely overabundant light WIMP with Q s td ^ 10 12 {m x / 
lOOGeV) 4 (if the production is thermal with chemical equilibrium as is 
usual). Underdense (ri st d < ^cdm) WIMPs have one or two solutions 
tt x = tt cdm per 77, if O std > l(r 5 (100GeV/m x ) and i) > 1(T 7 (100GeV/m x ) 2 



(^cdm/^std) (f° r 2rh > 5 MeV) 63]. In particular the neutralino density 



can be that of cold DM in almost any super symmetric model, provided 
10 12 (m x /100GeV) 4 > Q std > l(T 5 (100GeV/m x ) and the high energy the- 
ory accomodates the necessary combinations of values of bjm^ and Trh- 

Let us comment on other DM candidates. Sterile neutrinos u s would 
also be remnants of the pre-BBN era. If they are produced through oscil- 
lations with active neutrinos V g th eir pro duction rate has a sharp peak at 
T max ~ lSMeV^/leV) 1 / 3 ) [H, [3, M, Q which for m s > 1(T 3 eV is 
above 1 MeV. "Visible" v s (i.e. those that could be found soon in neutrino 
experiments) must necessarily have mixings sin(#) with v a large enough to 
be overabundant, and thus be rejected, in standard cosmologies. In LTR 
with Trh < Tmaxj the relic abundance of visible v s could be reduced enough 
for them to be cosmologically acceptable, both if they are lighter or heavier 
than lMeV 3^, 53, 67, 72]. E.g. for v s lighter than 1 MeV produced through 
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Fig. 7.6. The Hubble parameter H as a function of the photon temperature T 
before primordial nucleosynthesis for several cosmological models. 



53l . 1671 ] thus n s is small for 



oscillations, n s /n a ~ 10 sin 2 26 (T RH /5 MeV) 3 
low Trh, even if sin 6 is large. Another example is that of thermally produced 
axions, whose abundance can be strongly suppressed if Xrh is smaller than 
their freeze-out temperature ~ 50 MeV in standard cosmologies [j^ . 74]. 
Also super WIMPs may be produced in LTR models [7r3 |. 

Finally, let us remark that LTR scenarios are more complicated than the 
standard cosmology and no consistent all-encompassing scenario exists yet. 
In particular Baryogenesis should happen during the reheating epoch too, 
possibly through the Affleck-Dine mechanism 24j, |44j, 1471 . [52 ] . 



7.4-2 Models that only change the pre-BBN Hubble parameter 

We will consider two of these models, in which the change in WIMP relic 
density is more modest than in LTR: kination and scalar-tensor gravitational 
models. An homogeneous field (j), e.g. a candidate for quintessence, has an 
energy density = (j) 2 /2 + V ((/)). Kination is an epoch in which the kinetic 
term dominates over the potential V(<p) so ptotai — 4> 2 : /2 ~ a" 6 . No entropy 
is produced in this period, so T ~ a" 1 as usual. Thus H ~ ^/p t otai ~ ^ 3 
(see line "K" in Fig l7.6p . This case is intermediate between LTR, for which 
H ~ T 4 (see the line "LTR" in Fig |7.6|) and the standard radiation domina- 
tion case, for which H ~ T 3 (see the line "RD" in Fig 17. 6|) . Thus kination 
yields freeze-out temperatures T/. . larger than the standard, somewhere in 
between the LTR and the standard values. The only entropy dilution of the 
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density comes from the conversion of a larger number degrees of freedom 
present at the higher 7y o . into photon degrees of freedom at low tempera- 
tures, as particles annihilate and heat up the photon bath, and this effect is 
modest. The contribution of the <p kinetic energy to the total density is usu- 
ally quantified through the ratio of </> -to- photon energy density, i]^ = p^/ p 7 
at T ~ 1 MeV so that at higher temperatures H ~ y^(T/lMeV)i?staiidard- 
Notice that at T — 1 MeV, i.e. during BBN, the quintessence field cannot 
be dominant, thus r]^ < 1. Ref. [50] finds that the enhancement of the relic 
density of WIMPs in kination models is 

^kinationMtd ^ y^lO 3 (m x /100GeV) (7.34) 

Thus, WIMPs that are underdense in the standard cosmology could account 
for the whole of the dark matter. 

Scalar-tensor theories of gravity 30, 3^, 51, 71] incorporate a scalar field 



coupled only through the metric tensor to the matter fields. In many of 
these models the expansion of the Universe drives the scalar field towards 
a state where the theory is indistinguishable from General Relativity, but 
the effect of the scalar field changes the expansion rate of the Universe at 
earlier times, either increasing or decreasing it. Theories with a single matter 



sector typically predict an enhancement of H before BBN. In Ref. 5l[ the 
H is enhanced by a factor A, which is A ~ 2.19 x 10 14 (T /T) (T is the 
present temperature of the Universe) for large temperatures T > T^. At 
Tff,, A drops sharply to values close to 1 before BBN sets is (see the line 
"ST1" in Fig ESI). WIMPs freeze-out at T > T while H ~ T L2 , but at 
the transition temperature T^, H drops sharply to the standard value, and 
becomes smaller than the WIMP reaction rate. The already frozen WIMPs 
are still abundant enough at to start annihilating again. This is a post 
freeze-out "reannihilation phase" peculiar to these models. The WIMP relic 
abundance is reduced in this phase, but nonetheless remains much larger 
than in the standard case. The amount of increase in the WIMP relic 
abundance was found in Ref. 51] to be between 10 and 10 3 . With more than 
one matter sector, of which only one is "visible" and the other "hidden", 
scalar-tensor models may also produce a reduction of H by as much as 
0.05 of the standard value (see line "ST2" in Fig 17. 6p before the transition 
temperature at which H increases sharply to the standard value before 
BBN 71]. The maximum reduction of the WIMP relic abundance is larger 



for larger WIMP masses, ranging from a factor of 0.8-0.9 for masses close 



to 10 GeV to 0.1-0.2 for those close to 500 GeV [51] 
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